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ABSTRACT 

We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the 
X-ray binary Cygnus X-l in the hard state. This is the first time this state has been observed in Cyg X-l with 
NuSTAR , which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm 
that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires 
a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. 
We use the reflection models of |Garcia et al. ( 2014| ) to simultaneously measure the black hole spin, disk inner 
radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level 
of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner 
disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3a 
confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to 
hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather 
than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron 
line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection 
model. 
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1. INTRODUCTION 

The hard state in X-ray binaries (XRBs) is characterised by 
a spectrum with only weak thermal emission from the disk, a 
hard Comptonized component, and reflection from the accre¬ 
tion disk. It is generally agreed that the change between soft 
and hard states must be due to a change in the geometry of the 
accreting system, but there is still controversy over the exact 
nature of the change. This disagreement hinges on whether 
the accretion disk is significantly truncated in the hard state. 
It is generally agreed that the disk is truncated at a radius of 
^100 gravitational radii ( tq) as binaries go into quiescence 
(e.g. |Tomsick et al. 2009| ), but at what stage truncation be¬ 
gins as the accretion rate drops remains uncertain. In the disk 
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truncation model ( [Esin, McClintock & Narayan|1997[ ), as the 
accretion rate drops the classical thin disk becomes unstable 
and is truncated. This radius is typically greater than ~10- 
20 re, and the inner disk is replaced with a radiatively i neffi 
cient advection dominated accretion flow (ADAF; Narayan & 


1 1 994[ |Narayan| 1 99~6| |B landlord & Begelman|i999|). This 

haviour has been invoked to explain the transition from soft 


behaviour has been invoked to explain 1 
to hard states. However, some authors claim instead that the 
disk extends to small radii in the hard state while the accre¬ 
tion rate is higher than ~ 1 per cent of the Eddington rate, as 
measured using the profile of the relativistic iron line (e.g. 


(m3 

Re 


[Miller, Homan & Miniutti 


[Reis, Fabian & Miller|2010| . 


2006] IMiller et al.|[2006l [2015; 


There are two main ways of probing the inner edge of the 
thin disk, both of which rely on measuring the relativistic ef¬ 
fects that occur close to the black hole. The first involves 
modelling the relativistic blurring of the Fe Ka line emit¬ 
ted when coronal X-rays are reflected from the inner accre¬ 
tion disk ( [Fabian et al. fl 989| . This technique has been used 
extensively to measure black hole spins, b oth in XRBs and 
active galactic nuclei (AGN; see review by [Reynolds||2013 [ ). 
The second method, known as continuum fitting, relies on the 
black body spectrum of th e disk itself (see review by McClin 


tock, Narayan & Steiner||20l4]). There are difficulties with 
using either method to look at XRBs in the hard state, as both 
reflection and disk components are much weaker in this state 
than in the soft state. In this work we will focus on modelling 
the iron line, as the disk black body lies outside the energy 
range covered by NuSTAR. The expected behaviour of the re¬ 
flection spectrum, if the disk truncation model is correct, is for 
the iron K line to become narrower as the disk truncates, since 
the relativistic blurring becomes less extreme. While this is 
difficult to measure due to the relative weakness of the reflec¬ 
tion component in the hard state, truncation has been claimed 
in some sources at relatively high accretion rates (e.g. [Plant 
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|et al.12015] ). This also provides a natural explanation for the 
lower reflection fraction in the hard state, as far less of the 
Comptonized emission encounters the thin disk if it is trun¬ 
cated, assuming a compact corona. 

Further confusing the issue is the degeneracy between 
source hei ght and inner radius ( [Dauser et al. 2013 Fabian 


et al.||2014| ). For simple coronal models where the source is 


compact and on-axis (the ‘lamp-post’ geometry; see Wilkins 
|& Fabian|2012| for discussion of more complex geometries), 
the narrow lines predicted by truncation can instead be inter¬ 
preted as an increase in the distance between the source and 
the disk, which leads to less extreme light-bending and more 
illumination of the outer disk. 

Cygnus X-l (Cyg X-l) was the first confirmed black hole, 
when its mass w as determined to be too great for it to be a 
neutron star (e .g. Webster & Murdin 1972, Tananbaum et al. 


|1972[|Bolton||1972[ Gies & Bolton||1986| ). It was one of the 

first X-ray sources discovered, u sing Geiger counters aboard 
a sounding rocket (Bowyer et aLl |l965| ), and has been ex¬ 
tensively studied since then. Several authors have investi¬ 
gated the inner disk in the hard state in this source. |Done| 
& Zycki]( l999) used a combination of ASCA , EXOSAT and 


Ginga data to investigate the hard state, finding significant 
relativistic blurring, consistent with a moderately truncated 
disk. Young et_al. ( 2001| use data from the same three in¬ 
struments and argue in favo ur of a non-truncated d isk around 
a Schwarzchild black hole. Frontera et al.|( [2001| ) compared 
the soft and hard states using BeppoSAX data, finding inner 
radii of < lOrp for the s oft s tate and lOl^ro in the hard state. 


[Reis, Fabian & Miller] ( |2010| found an inner radius consistent 
with t he IS CO of a Schw arzchild black hole using Suzaku 
data. |Miller et al.| ( ]2012| ) found that the inner edge of the 
accretion disk in Cyg X-l was consistent with being located 
at the innermost stable circular orbit (ISCO) using reflection 
modelling of 20 Suzaku observations of the hard state, and 
|Fabian et al.[ ( |2012aj ) found that the relativistic reflection re¬ 
quired high spin, fitting an average observation from the same 
dataset. Measurements of the spin of Cyg X-l in the soft state 
using continuum fitting and and relativistic reflection have 
also r eturned high spin values ( |Gou et al.|201 l[|Tbmsick et ah 


2014), and the same has been found during transition (Duro 
et al.12011). Therefore, a systematic difference in the spin 


measurement between states would be indicative of trunca 
tion of the disk. 

Here we present NuSTAR and Suzaku observations of 
Cyg X-l in the hard state. The soft state of Cyg X-l was 


previously observed with NuSTAR and Suzaku (Tomsick et al. 
|2014| ), but this is the first time that the hard state of Cyg X-l 
has been seen with the ground-breaking sensitivity and band¬ 
pass of NuSTAR. The launch of the NuSTAR X-ray telescope 
(Harrison et al.|20131) has revolutionized the study of the X- 
XRBsfa 


ray spectra of! 


, as NuSTAR’s high sensitivity, broad en- 


precedented quality ( 

Miller et al. 

2013b a| Fiirst et al. 

2013 

|2014; Natalucci et a] 

. 2014 ; Tenc 

iulkar et al. 2014). 1 

n this 


paper we focus on modelling the broad band spectrum of the 
Cyg X-l hard state. Timing and variability analysis will ap¬ 
pear in future work. In particular, we are interested in mea¬ 
suring the spin and inner radius of the disk, so that we can es¬ 
tablish whether the disk is truncated, and if so by how much. 

2. OBSERVATIONS AND DATA REDUCTION 

Cyg X-l was observed for ^ 35 ks with NuSTAR and 
~ 107 ks with Suzaku. The details of all observations are 


Table 1 

Times, exposures and count rates for the two NuSTAR detectors (FPMA and 
FPMB) and three Suzaku detectors (XIS1, PIN and GSO). Exposure times 
are after filtering for background flares. Count rates are for the energy bands 
used for spectral fitting: 3-79 keV for NuSTAR ; 1.2-1.7 and 2.5-9 keV for 
XIS1; 20-38 and 43-70 keV for the PIN; and 60-300 keV for the GSO. 


Instrument 

Start Time 
(UT) 

On source 
exposure (ks) 

count rate 
(s- 1 ) 

FPMA 

2014 May 20th, 09hr 

34.4 

205 

FPMB 

” 

35.3 

189 

XIS1 

2014 May 20th, 07hr 

12 

134 

PIN 

2014 May 19th, 07hr 

107 

9.6 

GSO 

” 

107 

8.4 


given in Table [T] In this section we separately describe the 
observations ana data reduction for each instrument. 

2.1. NuSTAR 

The NuSTAR data (ObsID 30001011007) for the two fo¬ 
cal plane modules (FPMs) were reduced using the latest ver¬ 
sions of the NuSTAR data analysis software (NuSTARDAS 
vl.4.1) and CALDB (20140414). Cleaned event files were 
produced using the NUPIPELINE routine, and spectral prod¬ 
ucts with nuproducts. We used a 150" circular extraction 
region for the source spectrum, and a 100" region for the 
background, taken from an area of the detector not strongly 
contaminated by source counts (i.e. the furthest away part of 
the field from Cyg X-l). The source flux is at least an or¬ 
der of magnitude higher than that of the background over the 
whole band, so variations in the background with position or 
off-axis angle will not affect our results. We fit the NuSTAR 
data over the whole energy range of 3-79 keV. The NuSTAR 
spectra are binned to oversample the spectral resolution by a 
factor of 3, and to a minimum signal to noise ratio of 50 after 
background subtraction. At low energies in both detectors the 
oversampling restriction dominates (below ^ 20 keV), as the 
total number of counts is extremely large (the FPMA spec¬ 
trum alone contains almost 5 million counts). 

2.2. Suzaku 

For Suzaku , we reduce d the X-ray Imaging Spectrometer 
(XIS;|Koyama et al.|2007|) data from ObsID 409049010 using 
HE AS OFT v6.16 and the 2014 July 1 version of the calibra¬ 
tion files. The XIS coverage began on 2014 May 20, 6.7 h 
UT and ended on 2014 May 22 7.4 h UT. Although there are 
normally three operational XIS units, only XIS1 was on dur¬ 
ing this observation due to limitations on the available power 
during the late stages of the mission. We produced two sep¬ 
arate XIS1 event lists using aepipeline, one of the full 
observation and one only including the first half of the obser¬ 
vation (from the start time given above to 2014 May 21, 9.1 h 
UT), corresponding to the time of the NuSTAR coverage. We 
ran aeattcor2 and xiscoord to update the photon posi¬ 
tions. Due to the brightness of Cyg X-l XIS1 was operated 
in 1/4 window mode, which reduces pile-up at the expense of 
reduced exposure time, and used the 0.3 s burst option. Even 
with the window mode, we used pileest and found >4% 
pile-up within l f of the center of the point spread function. 
For making the XIS 1 spectrum, we used an annular extraction 
region with an inner radius of 1' and an outer radius of 4'. Af¬ 
ter this region is excluded, the photon averaged pile-up frac¬ 
tion is 1.6% and the pixel averaged fraction is 0.7%. A back¬ 
ground spectrum was taken from a rectangular region near the 
edge of the field of view. A portion of the source extraction 
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region includes inactive regions of the detector, and we ac¬ 
counted for this by adjusting the BACKSCAL keyword in the 
source spectral file. After creating the source and background 
spectral files, we used xisrmf gen and xissimarf gen to 
produce the response matrix. In the spectral fitting, we in¬ 
cluded the 1.2-1.7keV and 2.5-9.0keV ranges. There are 
known calibration uncertainties associated with the Si K-edge 
in the 1.7-2.1 keV band and we extend the excluded range up 
to 2.5 keV because of the presence of narrow features poten¬ 
tially due to calibration errors. We bin the spectra to oversam¬ 
ple the spectral resolution by a factor of 3, and to a minimum 
signal to noise ratio of 30 after background subtraction. The 
oversampling restriction dominates below ~ 8 keV, so the full 
instrumental resolution is maintained over the iron line band. 
We do not add any systematic errorto the XIS data. 

The hard X-ray detector (HXD; fIhkahashi et al.|20071 ) con¬ 
sists of Silicon PIN diodes covering an energy range of 12- 
70 keV and GSO scintillators covering an energy range of 
40-600 keV. Cyg X-l was observed at the XIS nominal po¬ 
sition with a net HXD exposure of 107.0 ks, and detected 
up to 400 keV by GSO. The HXD data was analyzed in 
the standard manner using the perl script hxdpinxbpi and 
hxdgsoxbpi for PIN and GSO respectively. The mod¬ 
eled non X-ray background (NXB) was taken from public ftp 
site^H and cosmic X-ray background (CXB) was also sub¬ 
tracted based on previous HEAO observations ( Gruber et al.| 
1999} for PIN. A 1 per cent systematic error was added for 
eacn spectral bin. In the spectral fitting, we used the energy 
responses: ae_hxd_pinxinomell_20110601.rsp for PIN and 
ae_hxd_gsoxinom_20100524.rsp with an additional correction 
file (ae_hxd_gsoxinom_crab_20100526.arf) for GSO. We fit 
the PIN data from 20-70 keV, excluding lower energies to 
avoid thermal noise and the 38-43 keV band to avoid a known 
calibration feature, and the GSO data from 60-300 keV. We 
extracted spectra for both HXD detectors corresponding to the 
intervals simultaneous and non-simultaneous with the NuS¬ 
TAR exposure. We find no significant differences in the spec¬ 
tral shape for either detector between the two intervals, so we 
use the full exposures in all fits. 


3. RESULTS 

We split the spectral fitting into three sections: in we 
fit simple continuum models to the broad-band spectrum to 
establish the presence and energy of the high-energy cut-off 
without modelling the reflection spectrum; in §3.2| we fit the 
iron line profile in detail, to measure the relativistic blurring 
without c omp lications from other regions of the spectrum; fi¬ 
nally, in §3.3| we fit the full broad-band spectrum, including 
the reflection, to calculate our final best fit model. We use 
Xspec ( |Amaud|1996| ) version 12.8.2a for all fitting. All er¬ 
rors are quoted at 1 standard deviation unless otherwise stated. 


3.1. Broad-band continuum fitting 

To establish a baseline continuum model from which to 
measure the reflection properties we initially fit the NuSTAR 
FPMA and FPMB and Suzaku XIS, PIN and GSO data si¬ 
multaneously. For this section only, we include a 5 per cent 
systematic error to each spectral bin for all datasets, to en¬ 
sure that the model is being fit to the broad-band spectrum 
rather than being dominated by small but high signal-to-noise 

15 ftp://legacy.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.2_tuned/ 
and ftp ://legacy.gsfc.nasa.gov/suzaku/data/background/gsonxb_ver2.6/ 


Table 2 

Fit parameters for the continuum models (see text). We include 5 per cent 
systematic error on all continuum only fits, so x 2 values do not indicate 
absolute goodness of the fits. Cross-calibration constants are given for the 
best-fit model. 


Parameter 

Value 

Unit/description 

POWERLAW 

r 

1.702 db 0.003 

Photon index 

Tiisk 

0.43 ±0.01 

disk temperature (keV) 

X 2 /d.o.f. 

1973/948 


CUTOFFPL 

r 

1.568 ± 0.005 

Photon index 

E'cut 

156 ±3 

High energy cutoff (keV) 

Tiisk 

0.53 ±0.01 

disk temperature (keV) 

X 2 /d.o.f. 

941/947 


COMPTT 

^plasma 

43.4 ±0.8 

Plasma temperature (keV) 

Tlisk 

0.448 ± 0.007 

disk temperature (keV) 

T 

1.27 ±0.03 

Optical depth 

X 2 /d.o.f. 

863/947 


Cross-normalization constants 


Cfpma 

1 


Cfpmb 

1.015 ±0.004 


Cxisi 

0.948 ± 0.006 


CpiN 

1.31 ±0.01 


Cgso 

1.16 ±0.02 



features in the NuSTAR spectra. We stress that the x 2 val- 
ues presented in this section should therefore not be taken as 
absolute indications of the fit quality, only relative - without 
the systematic error, the reduced x 2 values are on the order 
of 35-40. In all fits we include interstellar absorption, mod¬ 


elled with TB abs (using the abundances of Wilm s, Allen & 
McCray 2000), fixed at the value of 6 x 10 21 


cm " from Tom 
sick et al" ( |2014| . We u se o nly t he simultaneous XIS1 for the 
broad-band fits (§ |3.1| and |3.3|), as this ensures a consistent 


power law index T between the XIS and FPM spectra. The 
full XIS1 spectrum gives a significantly different value of T, 
making it much more difficult to simultaneously fit the XIS 
and NuSTAR spectra. We include a constant normalization 
offset between the instruments, to allow for differences in flux 
calibration. 

In the left panel of Fig. [I] we show the five spec¬ 
tra fit with a Comptonization plus disk black body model 
(tbabs*(diskbb+COMPTT)), along with background spec¬ 
tra and residuals. Clear reflection features are apparent in the 
spectrum, along with the absorption line at ^ 6.7 keV previ¬ 
ously seen by Tomsick et al. ( 2014| from the companion stel¬ 


lar wind. In thi"right panel", we show the unfolded spectrum 
to this model, where it becomes immediately apparent that 
there is a steep turnover at high energies, just above the NuS- 


with RXTE and Suzaku (e.g. 

Gilfanov, Churazov & Revnivt- 

sev 2000 Wilms et al. 2006; 

Nowak et al. 201 lj>, and appears 

to be similar between observations. The presence of this high 


energy curvature, while not immediately apparent in the NuS¬ 
TAR data alone, must be taken into account befo re the reflec¬ 
tion spectrum can be accurately modelled (§ |3.3| ). The model 
parameters for this fit are give in Table [2j along with those 
for simpler power law models (with anawithout a cut-off). 
We find the best fit with the COMPTT model, which we use as 
o ur in itial continuum model for the detailed spectral fitting in 

m 

We also fit the more complex EQPAIR Comtponization 
model to these data, however with the fitting procedure out- 
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Figure 1 . Left : Source and background (shaded regions) spectra from the five detectors, fit with a Comptonization model (COMPTT; |Titarchuk|1994|, and 
residuals to the fit. A constant offset is allowed between the data sets to account for flux calibration differences. Right: The same five spectra, shown unfolded 
against the model, showing the sharp high energy turn over. XIS1 data are rebinned slightly in XSPEC for clarity. 


lined above this results in unphysical best-fit parameters, as 
the model is complex enough to be partly able to fit the reflec¬ 
tion curvature as well as the continuum. We therefore defer 


discussion of this model to the more detailed fitting in § 3.3 


The accuracy and reliability of the cut off measurement 
depends on the accuracy of the GSO background. From 
~ 100 keV upwards the GSO spectrum is background dom¬ 
inated, so small differences between the modelled and true 
background can potentially have a large effect. The back¬ 
ground spectrum itself is completely dominated by the non 
X-ray background (NXB), genera ted by radioactive iso topes, 
cosmic rays, and albedo neutrons ( [Kokubun et al.|2007] ). This 
background signal is significantly time variable, however it 
is well understood and can be accurately modelled and com¬ 
pared with the observed spectrum during Earth occultations. 
The estimated systematic error on the background spectrum is 
less than 1% for observations of over 10 ks ( |Fukazawa et al.| 
2009), which we exceed by a factor of 10. To test the sys¬ 


tematic error this uncertainty introduces, we assume a con 
servative 1% error in the normalization of the GSO back¬ 
ground, calculate background spectra scaled up or down by 
this amount, and then re-fit the data using these background 
spectra instead. This changes the cut off energy by ±1.5 keV 
(comparable to the ±0.8 keV statistical error), but makes no 
qualitative difference to the spectral shape. 


3.2. Iron line fitting 

To avoid possible confusing effects from other spectral fea¬ 
tures we fit the XIS and NuSTAR spectra over the iron line 
band (4-10 keV) alone. For this section we use the full XIS 
exposure, rather than just the strictly simultaneous data, as 
changes in the continuum shape are less critical when fitting 
a limited energy range, and we want to maximize the signal 
around the iron line. There is a slight difference in the photon 
index measured by Suzaku and NuSTAR when using the full 
dataset, which is due to spectral variability (the simultaneous 
Suzaku and NuSTAR data agree well, and the two Suzaku spec¬ 
tra from simultaneous and non-simultaneous intervals show 


the same difference). We therefore allow the index to change 
between the two, though this is a small effect (~ 2 per cent), 
and will not affect our results or conclusions. 

In Fig. [2] we show the iron line profile from 4-10 keV, fit 
with several different models. In the top panel we show the 
total profile, fit with a powerlaw (including interstellar ab¬ 
sorption). The line profile is clearly broad, extending from 
~ 5-8 keV, with an absorption feature at around 6.7 keV. 
There is a generally very good agreement between the XIS 
and FPM spectra, although there is a difference in flux be¬ 
tween the two above ~ 8 keV, which holds whether we use 
the full XIS spectrum or just the simultaneous data, implying 
it is due to calibration differences. There are known calibra¬ 
tion issues with XIS1 spectra above ^ 8 keV, stemming from 
variable background contamination from cosmic-ray photons 
( Ishida et~aLl|2011| ), making the background level uncertain. 
As NuSTAR is likely to have better calibration at high ener¬ 
gies and dominates the signal in the combined spectrum, this 
disagreement should not significantly affect our results. We 
next fit a series of models to the line profile, shown in the re¬ 
maining panels of Fig.[2| and the model parameters are shown 
in Table [3] Before fitting a blurred reflection spectrum, we 
test simpler models to check whether the line profile can be 
explained as a combination of narrow features and absorp¬ 
tion, as the Cyg X-1 binary system is known to have a strong 


stellar wind which can produce these effects (Marshall et al. 


2001; Schulz et al. 2002] |Miller et al.|2005| ). We first fit a sim- 
ple model consisting of a powerlaw plus two narrow emission 
lines (POWERLAW+GAUSS+GAUSS), but this model leaves a 
strong, broad excess from 5-6 keV. Modelling the line with 
ionized, but not blurred, reflection (using XlLLVEE 


16 


Garcia 


et al. 2013| ) fares no better, even when we include absorption 


from the wind using the same xstar (Kallman & Bautista 


16 We note that Xillver version 0.2d contains a bug affecting models 
with cut-off energies > 300 keV, which affects both Xillver and Relxill 
(Garcia & Dauser, private communication). This is not a problem for this 
work, as all of our models have cut-off energies < 300 keV. 
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Figure 2. Data to model ratios for the FPMA, FPMB 4-10 keV spectra and 
XIS1 4-9 keV spectrum. In all cases we include interstellar absorption. Num¬ 
bers in brackets correspond to model numbers in Table[3] 


|2001| ) grid as |Tomsick et al.| ( 2014| ). 

Tne exceptional quality or the combined spectrum makes 
it clear that only a relativistically blurred iron line can ex¬ 
plain the observed line profile. The broad excess remain¬ 
ing from 5-6 keV after fitting the two Gaussian model must 
be due to blurred reflection, as there are no strong atomic 
lines in this band. We fit the_relativistic reflection spectrum 
with the relxilllp model ( [Garcia et al.||2014| ), which self- 
consistently calculates the reflected emission as a function of 


Table 3 

Fit parameters for the model fits to the iron line band (4-10 keV). Model 1 is 
a pure power law model, model 2 includes two narrow Gaussian lines, 
model 3 replaces the lines with the xillver reflection model, model 4 
includes an ionized absorber, modelled with an XSTAR grid, and model 5 is 
the full relativistic reflection model (see text for description). All models 
include interstellar absorption. Errors are only quoted when x? < 2. We 
give the cross-calibration constants for the best-fit model 


Model 

Parameter 

Value 

Units 

1 

r 

1.85 



X 2 /d.o.f. 

4311/256 


2 

r 

1.85 




6.34 

keV 


Wx 

33 

eV 


E 2 

7.06 

keV 


w 2 

20 

eV 


X 2 /d.o.f. 

1574/251 


3 

r 

1.87 



Ape 

1.01 



logCref 

e 

X 2 /d.o.f. 

1.67 

68.1 

1694/251 

log(erg cm s -1 ) 
degrees 

4 

Nxx 

6.3 x 10 21 

cm -2 


l°g 'Cabs 

3.29 

log(erg cm s _1 ) 


r 

1.88 



^Fe 

1.01 



log Cref 

1.67 

log(erg cm s _1 ) 


e 

X 2 /d.o.f. 

68.1 

1395/249 

degrees 

5 

nxx 

log Cabs 

1A7 -0A7 X 10 22 
4.85 ±0.02 

cm -2 

log(erg cm s -1 ) 


h 

< 1.21 

Hsco 


a 

> 0.84 



n n 
r 

Ape 

1.11 ±0.01 

1 926 +0 ' 001 

±.»^u_ 0 006 

< 0.55 

nsco 


logCref 

2.768 ± 0.003 

log(erg cm s _1 ) 


0 

47.8 ±0.3 

degrees 


R 

0.809 ±0.01 



-F/Gauss 

6.43 ±0.01 

keV 


FkGauss 

12 ± 1 

eV 


X 2 /d.o.f. 

294/249 


Cross-normalization constants 



Gfpma 

1 



Cfpmb 

1.0157 ±0.0007 



Gxisi 

1.148 ±0.008 



inclination for an on-axis point source above the disk. Rela¬ 
tive to models which assume a broken power law emissivity 
profile, this has the advantage of parametrizing the profile in 
terms of a physical quantity, reducing the number of parame¬ 
ters needed to describe the profile from 3 to 1, and restricting 
the profile to regions of parameter space that make physical 
sense. In addition to the broad line, we find a weak narrow 
emission feature at 6.4 keV which we model with a narrow 
Gaussian line (see bottom 2 panels). This feature is only ob¬ 
vious in the full XIS spectrum, but corresponds to an excess in 
the NuSTAR spectrum as well, and results in a highly signifi¬ 
cant improvement in the fit statistic A% 2 = 56 for 2 additional 
degrees of freedom. The iron line in Cyg X-l was shown to 
be a composite of broad and narrow lines using the Chandra 
Hi gh Energy Tr ansmission Grating Spectrometer (HETGS) 
by Miller et al . d2002| . The best fit model parameters are 
showiTin Table 0 

In Fig. [ 3 ] we show A% 2 contours from fitting the iron line 
profile for the spin (a), source height ( h ), and inner radius 
(tin) model parameters. All three are strongly constrained to 
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Figure 3. Left, x 2 contour plot for the inner radius of the accretion disk, in units of the ISCO. Middle : % 2 contour plot for black hole spin parameter. Right: 
contour plot for the height of the emitting source above the disk plane, in units of the ISCO radius. All plots are calculated using the STEPPAR function in Xspec. 


values indicating strongly blurred emission from the inner¬ 
most regions of the accretion disk. Even if we assume the 3a 
limits on a and tin (0.8 and 2 nsco) this gives a inner radius 
of only ~ 6ro, which argues strongly against significant disk 
truncation. 

We check for degeneracy between parameters by using a 
Markov chain Monte Carlo (MCMC) algorithm to explore the 
parameter spaced] We use 50 walkers with 15,000 iterations 
each, burning the first 5,000. The walkers are initially dis¬ 
tributed around the best fit value found by Xspec, following 
a normal distribution in each parameter, with the standard de¬ 
viation set to the delta value of that parameter (we ensure that 
the delta values are all smaller than the uncertainties on each 
parameter). We check that the auto-correlation time is at least 
100 times smaller than the chain length for all parameters, and 
‘jack-knife’ the chains (comparing the distributions from the 
first and second halves of the chain) to check that there are no 
large systematic differences between the distributions. 

We calculate the distributions of all parameters (Fig. |4f^| 
and find significant degeneracies in only two cases. The trun¬ 
cation radius of the disk is degenerate with the spin, as ex¬ 
pected as both parameters control the effective inner radius. 
While these parameters cannot therefore be independently 
constrained, this remains a very strong limit on the truncation 
of the inner disk. If we take the 3a limits of the distribution, 
at no point does the inner radius exceed 3ro, so this strongly 
requires the disk not to be truncatecp) 

Additionally, the reflection fraction and iron abundance are 
degenerate. This is also expected since we are only fitting the 
iron line band, so a true constraint on either value is hard to 
achieve - within this limited energy range, the main effect of 
both of these parameters is to alter the strength of the iron 
line. This goes some way toward explaining the sub-solar 
iron abundance we find, which is unlikely to be correct as the 
iron abundance in Cyg-X-1 is generally thought to be super¬ 
solar (|Hanke et al.|2009[|Duro et al.|20 1 1 [ Fabi an et al.|2012bt 


Tomsick et al. 2014). We note that there is a slight discrepancy 
betweenthe steppar and MCMC contours - for example, the 
3-a limit on a is lower when calculated using STEPPAR. This 
difference can be explained by the limited accuracy possible 
with STEPPAR - while both methods agree on the requirement 


17 We use the xspec. em cee code by Jeremy Sa nders, based on the em¬ 
cee python implementation (Foreman-Mackey et al. 2013) of the Goodman- 
Weare affine invariant MCMC ensemble sampler (Goodman & Weare 2010). 

18 This plot and the plot in Fig. [ 7 ] are prod uced using the Triangle.py 
python module (Foreman-Mackey et al.|2014| 

19 This differs from the 6 rQ previously quoted because that estimate did 
not take into account parameter degeneracies 


for high spin, the best fit found by STEPPAR has a % 2 value of 
298, compared to the solution with 294 found by the MCMC 
analysis. In addition, the power law index and ionization pa¬ 
rameter £ re f are degenerate, although this covers only a very 
small range in log £ re f. 

To investigate the iron abundance further, we re-fit the data 
with the abundance fixed at 2, as found by Tomsick et al 

/oni Jk Uic _ 1 C anr 


(2014). This worsens the fit significantly (A%~= 15 for one 
fewer degree of freedom), however it introduces only minor 
differences in the other parameters (all of the relativistic blur¬ 
ring parameters are consistent within errors). The only param¬ 
eter to change significantly is the reflection fraction, which 
decreases from 0.81 db 0.01 to 0.34 =b 0.01. Because of the 
lack of change in other parameters and absence of significant 
residual features it seems likely that the constraint on the iron 
abundance is being driven by random statistical fluctuations 
in the data or a minor feature of the reflection model. 

In the top panel of Fig. [5] we show a comparison of the iron 
line profiles from the soft and hard states. We use the soft state 
observations 30001011002 and 30001011003, which were 
presented by Toms ick et al. ( 2014|), and were ta ken at a simi¬ 
lar orbital phase (0.845-0.962 in Tomsick et al. compared to 
0.902-0.104 in this work, so the absorption should be simi¬ 
lar). We plot the combined FPMA and FPMB spectrum from 
both observations. These data are reduced and extracted as in 
Tomsick et al. (2014] ). We fit both spectra with simple contin¬ 


uum models over the 3-4 and 8-10 keV bands. For the hard 
state, we use a power law, and in the soft state a power law 
plus disk b lack body, with the tem perature fixed at the best fit 
value from Tomsic k et al. | ( |2014] ) . We then scaled the residu¬ 
als of the hard state by a factor of 3 so that the amplitudes of 
the two lines are comparable. It is immediately apparent from 
the residuals that there are no significant qualitative changes 
to the red wing of the line, which is the part of the line profile 
which should be most affected by truncation of the disk. In 
both cases, the line extends to below ~ 5 keV. The ~ 6.7 keV 
absorption feature appears to be consistent between the two, 
with the spectral changes concentrated on the blue side of the 
line. In the lower panel we show the same profiles, but after 
correcting for the absorption from the companion wind (us¬ 
ing the absorption parameters from model 5). This makes the 
6.4 keV excess from the narrow line in the hard state more 
obvious, and demonstrates more clearly the lower energy of 
the peak in the hard state. 


3.3. Full spectrum fits 

We now fit all five spectra (XIS, FPMA, FPMB, PIN and 
GSO) simultaneously over the full energy range, modelling 
both the continuum and reflection components. We allow for 
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Figure 4. Output distributions from the MCMC analysis of the best fit model of the iron line profile. Contours correspond to 1, 2 and 3 a. Points are shown 
where the density of points drops beyond the 3cr limit, y-axes for the histograms are in arbitrary units. Only two of the contour plots show strong degeneracy: 
a and r m are degenerate, as are the iron abundance and reflection fraction. Negative values for ri n and h correspond to units of nsco> following the convention 
from Relxill. 


a constant normalization offset to accoun t for flux calibration 
differences between the spectra. As in § |3.1j we use only the 
XIS1 data that was taken simultaneously with NuSTAR. 

We initially fit the full spectrum with a combinati on o f the 
models for the continuum (COM PTT + DISKBB; § ED and 
line profile (Relxill + gauss; §|3.2|). The full model is, in¬ 


cluding absorption, TBABS * XSTAR* (DISKBB + COMPTT + 
relxillp + GAUSS). As the relxill models assume a cut¬ 
off power law as the illuminating spectrum, we fix the cut-off 
of the reflection spectrum to three times the plasma temper¬ 
ature of the Comptonization component (i.e. the peak of the 
Wien law in energy flux), which is approximately the cut-off 
energy that would be measured with a power law model. This 
model does not result in a formally acceptable fit - we find no 
solutions with xl < 2. It is not clear why this model does not 
work, as it results in significant residuals over the full energy 
band. We conclude that it is likely due to subtle spectral cur¬ 
vature that cannot be accounted for by the simple continuum 


model used. We test adding a partial-covering ionized absorp¬ 
tion zone as a potential explanation for the curvature, and find 
that while it significantly improves the fit (A % 2 ^ 500, for 
3 additional degrees of freedom) the fit is still not acceptable 
(xl = 1669/839 = 1.99) 

Next, we fit the more sophisticated EQPAIR model ( |Coppi[ 
1999) . EQPAIR calculates the Comptonized spectrum for a 


combined population of thermal and non-thermal electrons. 
EQPAIR was used extensively to_model the broad band spec 
trum of Cyg X-l by |Nowak et al. (2011), who found that it 


could account for the spectral"shape well. As in Nowak et al. 
we assume the default distribution of non-thermal electrons 
a density following a power law 7“ 2 between 7 m i n = 1.3 and 
7 max = 1000, where 7 is the Lorentz factor. While this model 
can also calculate the reflection spectrum as of version 1 . 10 , 
we set the flux of this component to zero and continue to use 
RELXILLLP for our reflection model. This is done for two rea¬ 
sons: firstly to maintain consistency with the other models and 
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* l.i- 


Energy (keV) 


Figure 5. Top : Comparison of the NuSTAR iron line profiles of Cyg X-l in 
the hard and soft states, with observations taken at approximately the same 
orbital phase. The hard state profile is scaled to approximately match the flux 
of the soft state line at 6 keV. Bottom-. As top, but after correcting for the 
^ 6.7 k eV absorption line using the parameters from the best fit models in 
jTomsick et al. 1^2014) and this work. 


the fits in § 3.2 and secondly because the relxill models 
are more sophisticated and include more atomic physics than 
the ireflct model used by eqpair. We fix the cut-off of the 
illuminating spectrum of the reflection model to 100 keV, as 
there is no parameter to tie it to in the EQPAIR model, and it is 
very hard to constrain independently. We choose 100 keV as 
it corresponds to the peak of the spectrum shown in Fig.H] and 
is in the range of 2-3 times the plasma temperature (Table [2]) 
expected for a powerlaw approximation to a Comptonization 
spectrum. While this is not an exact value, the systematic 
error introduced is unlikely to be large, comparable to that in¬ 
troduced by the reflection model assumin g a po wer law as the 
input spectrum. Following [Nowak et al.| ( |201 f) , we include an 
un-Comptonized black body as well as the Comptonized one 
included in the eqpair model. This model can better account 
for the subtle spectral curvature, and results in a much bet¬ 
ter fit to the broad band spectrum (x 2 = 1219/919 = 1.33). 
The data, model and residuals are shown in Fig. [6j and the 
parameters are given in Table [4] 

We note that in the best fit model almost all the black body 
flux is outside the range covered by the regions of the XIS1 
spectrum that we fit. It is therefore likely that the tempera¬ 
ture is being affected more strongly by the EQPAIR contin¬ 
uum component than the black body itself. The measured 
temperature, while consistent with the hard state in Cyg X-l 
( [Balucinska-Church et al.|1995||Di Salvo et al.|2001| ), should 
therefore be treated with a degree of caution. 

We show the confi denc e contours for the reflection spec¬ 
trum in Fig. [7] As in jj |3.2| these are calculated using MCMC 
chains. As we are exploring a larger parameter space, we in¬ 
crease the number of walkers to 100, and the chain length to 
20,000. The results are very similar to those found from fitting 


Table 4 

Fit parameters and normalization constants for the best fit broad band 
model. Zh2s, and l nt refer to the hard, soft, and non-thermal compactnesses, 
respectively. The X-ray luminosity is absorption corrected, however the 
fluxes are not. We assume the distance of 1.86 kpc from Reid et al. ( 2011 > 
and mass of 14.8 Mq from Orosz et al. ( 2011). The reflection traction R is 
calculated from the ratio of model fluxes, rather than the model parameter, 
hence no error is calculated. 


Parameter 

Value 

Units 

Interstellar absorption (tbabs) 

Nh (5.1 ± 0.1) x 10 21 

cm -2 

Wind absorption (xstar) 
Nh 

(1.35 ± 0.07) x 10 22 

cm -2 

k>g(£abs) 

4.81 ±0.02 

log(erg cm s 1 ) 

disk black body (diskbb) 
±bb 

O 

O 

1 + 

O o 

o o 

O l- 1 

keV 

Comptonization (eqpair) 

k /Is 
lnt / 

0 7^+ 0 - 06 

U -' °-0.08 

0 O60 +0 ' 003 
u.uou_ 0 010 


r 

0.46 ±0.01 


Narrow line (GAUSS) 

-^Gauss 

6.40 ±0.01 

keV 

H^Gauss 

10 ± 1 

eV 

Reflection (relxilllp) 
h 

< 1.56 

nsco 

a 

> 0.97 


An 

1.5 ±0.3 

nsco 

<9 

45.3 ±0.4 

degrees 

r 

1.60 ±0.02 


k>g(£ref) 

^Fe 

3.03 ±0.01 

4.7 ±0.1 

log(erg cm s _1 ) 

R 

0.16 


Cross-normalization constants 


C'fpma 

1 


Cfpmb 

1.0163 ±0.0005 


C'xisi 

0.923 ± 0.002 


Cpin 

1.295 ±0.003 


C'gso 

1.192 ±0.005 


Fluxes 

2-10 keV (XIS1) 

7.57 x 10“ 9 

erg cm -2 s -1 

2-10 keV (FPMA) 

6.35 x 10“ 9 

erg cm -2 s -1 

10-50 keV (FPMA) 

1.04 x 10“ 8 

erg cm -2 s -1 

Tx, 0 . 1 - 200 /^Edd (FPMA) 

0.0250 



the line profile alone - the blurring parameters are strongly 
constrained to produce a high level of blurring, with a small 
inner radius, low source height, and high spin. Again, the 
spin and inner radius are highly degenerate, but are inconsis¬ 
tent with a large degree of disk truncation, with r- m < 3 vq in 
all cases. 

The EQPAIR continuum model includes a Comptonization 
of a population of non-thermal electrons. We test whether a 
hybrid electron population is important to the fit by fixing the 
compactness of non thermal electrons, / nt , to zero. We find 
a significantly worse fit with this model (Ax 2 = 222, for 
one additional degree of freedom), which indicates that non- 
thermal scattering is important in the continuum modelling. 


4. DISCUSSION 


As found by Tomsick et al. (2014), the inclination values we 
measure based on the reflection spectrum are not consistent 


with the bi nary inclination of 27.1 ± 0.8 degrees from Orosz 
et al.| ( [20TT] ), based on optical measurements. |Tomsick et al. 


use a selection of models and find inclinations in the range 


42-69 degrees. Progressing from the recent work of Tomsick 


et al., we are using the latest reflection models which self- 


consistently calculate the reflection spectrum and relativistic 
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Figure 6 . Left : Upper panel shows the broad-band spectrum and best fit EQPAIR model; lower panel shows the ratio of the data to the model. Shaded regions 
show the background spectra for each instrument. Right : The best fit model, showing the different model components. 


blurring as a function of angle ( Garcia et aL| |2014 ), but we 
nevertheless find a similarly high inclination of 45.3 =b 0.4 
degrees from our best fit broadband model. Tomsick et al. 
discuss the high ionization value found in the soft state as a 
possible cause of higher inclinations, due to additional Comp¬ 
ton broadening of the line not accounted for by the model. 
This is unlikely to be the same in the hard state, which we 
find has ~ 20 times lower ionization than the soft state, so 
we would expect significantly less broadening and a corre¬ 
spondingly lower inclination value. Including the expected 
broadening 


a 

E 


V rn Q c z J 


( 1 ) 


for a given disk tem perature Tbb and energy E (Pozdnyakov, 
|Sobol & Syunyaev|1983[ ), does not significantly affect either 


the quality of the fit, or the inclination. As a high inclina 
tion value has now been found with NuSTAR in two observa 
tions with different reflection models and in different states 
this rules out a range of systematic effects, making it more 
likely that there is a real physical effect behind this measure¬ 
ment. The recedi ng jet in Cyg X-l has not been detected ( Stir¬ 
ling et al.] |2001| ), which points to a small inclination angleT 


however the jet angle, which would offer an independent con¬ 
straint, is not strongly constrained. It is possible that the dif¬ 
ferences between the binary inclination and our inclination are 
due to precession of the disk or a warp in the disk. A strong 
modulation with a peri od of 142 zb 7 days is id entified as the 
effect of precession by |Brocksopp et al. |1999 ), with an angle 
of up to 37 degrees between the binary plane and the disk. 
This would be easily enough to explain the observed discrep¬ 
ancy, however there are uncertainties in the precession mea¬ 
surement (jZdziarski, Pooley & Skinner||2011|), and if this is 
the case we would expect the measured inclination to change 
with this period, which has not so far been observed. Addi¬ 
tionally, this period has only been observed in the hard state 
so may not be able to explain the soft state result. A potential 
problem with a warped disk is that thick disks are not expected 


to align (Ivanov & Illarionov| 1997t[Fragilo|2009| ), so thi 
not be possible, depending on the state of the accretion 


this may 
flow. 

We find some minor differences in the results when we fit 
the Fe Ka line alone compared to when we fit the full spec¬ 
trum. The T value found for the reflection component in the 
br oad band model, where it is independent of the continuum 
(§ |3.3| ), is significantly lower than that fo und when fitting a 
power law to the low energy continuum (§ |3.2| ), although this 
does not appear to impact the other reflecti on p arameters. If, 
as has been suggested elsewhere ( Fabian et al.| 2Q14[ ), the con¬ 
tinuum spectrum seen by the disk differs from that seen by 
an external observer, then tying the T values of the reflec¬ 
tion and continuum components may not be appropriate (this 
has also been observed in GX 339-4 by Fiirst et al., submit¬ 
ted). Additionally, we find a large difference in the iron abun¬ 
dance between the two fits, with the line fit returning a very 
low abundance (< 0.51 solar) and the broad band fit giving a 
super-solar value (4.7 db 0.1). The estimate from the line pro¬ 
file fit cannot be treated as reliable, as a measurement of the 
iron abundance requires additional reflection features to estab¬ 
lish the relative strength of the iron line. The result from the 
broad-band fit is more plausible, but could still be affected by 
systematic effects such as the assumed density of the dislpq 
We note additionally that neither of the iron abundanc es fo und 
in this work is consistent with the value from ITomsick et al.l 
( ]2014| ). This difference cannot be due to changes in the iron 
abundance of Cyg X-l, and must therefore be due to system¬ 
atic effects unaccounted for by the modelling. 

The strength of the relativistic blurring in the hard state 
spectrum argues strongly against the disk being truncated. 
The line profile is not significantly narrower in the hard state 
compared to the soft state, and our best fitting models exclude 
an inner disk radius larger than 3ro at 3cr, which does not al¬ 
low for significant truncation in the switch from soft to hard 

20 Changes in the density can potentially alter the strength of the iron line 
relative to the soft excess and Compton hump without otherwise changing the 
shape of the line profile. A discussion of this and other density effects will be 
presented in future work (Fabian et al., in prep.) 
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Figure 7. Output distributions from the MCMC analysis of the best fit model of the broad band spectrum of Cyg X-l. Contours correspond to 1, 2 and 3 a. 
Points are shown where the density of points drops beyond the 3cr limit, y-axes for the histograms are in arbitrary units. Negative values for r- m and h correspond 
to units of nsco, following the convention from Relxill. 


states. This holds whether we fit the whole spectrum, or re¬ 
strict our fits to the line profile. A similarly constant line 
profile was found by |Reis et al.l (|201 1 ) in XTE J1752-223, 
where the shape of the line profile was not significantly af¬ 
fected by the state transition. The main difference between 
the hard and soft line profiles appears to be attributable to 
changes in the ionization parameter of the underlying reflec¬ 
tion spectrum, rather than the extent of the relativistic blur¬ 
ring. The blurring parameters we find in this work are gen¬ 
erally consistent with those found by Tomsick et al. ( 2014| , 
but the ionization parameter differs by an order of magnitude 
and can strongly affect the shape of the line p rofile, extend¬ 
ing it to higher energies (e.g. |Ross & Fabian||2005| ). If the 
disk is not truncated, another explanation for the lower reflec¬ 
tion fraction in the hard state is needed: the ratio of reflected 
to continuum fluxes in the 20-40 keV band is 0.16 from our 
best fit ting model, and 0.52 for model 7 from To msick et al.| 
(2014 )). Potential explanations include outflowing material 
(Beloboro dov||1999[ |Miller et al.||2012| ) or increased Comp- 


tonization of the reflected emission by the corona ([Wilkins & 
Gallo |2015| ). We find some changes in the shape of the iron 


line profile between the soft and hard states, although we do 
not make a detailed comparison. These differences are largely 
on the higher energy side of the line - the red wings of the two 
line profiles overlap very closely when the hard state residuals 
are scaled to match those of the soft state. A more detailed in¬ 
vestigation of the changes between the states will be presented 
in in future work. 

Measuring the spin or inner radius of XRBs in the hard state 
is significantly harder than in the soft state due to the reduced 
flux and lower reflection fraction. For this reason there are 
relatively few estimates of these parameters to compare our 
results to, and no strong consensus on true values. Our result 
(a > 0.97) is consistent with recent measurements of the spin 
in the hard state by Miller et akj ( |2012| ) (0.6 < a < 0.99) and 
Fabian et al. (2012b) (0.97 ^q.o2 )• Previous works on the hard 


state have been unable to simultaneously measure the spin and 
inner radius, instead using one as a proxy for the other, due to 
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the strong degeneracies and limited spectral quality available. 
Our new measurement thus represents a significant step for¬ 
ward in understanding the behaviour of Cyg X-l in the hard 
state, made possible by high quality NuSTAR spectra and the 
latest generation of reflection models. These spin values are 
also consistent with the near maximal values found in the soft 
state, where the line profile is clearer and the spin can be inde¬ 


pendently measured using the continuum fitting method (Gou 
let al.|201U|20l4l|Tomsick et al.|2014| ). 

We use the relxilllp relativistic reflection model 
throughout this work. Unlike other such models, RELXIL¬ 
LLP explicitly assumes a particular geometry - that of an on- 
axis point source. In reality, the corona with have some radial 
and vertical extent, so eventually this assumption will break 
down. However, this should not affect the results presented 
here, as the strength of the relativistic blurring is very high. 
The combination of high spin, low source height and small 
inner disk radius all point toward a compact source, close to 
the event horizon. Under these circumstances the point source 
assumption is valid beyond the limits of current instrumenta¬ 
tion. Introducing either a large radial or vertical extentto the 
source would flatten the em issivity profile ( Wilkins & Fabian 


2Q12| pauser et al.||2013) , reduce the level of blurring and 


worsening the fit. One possibility for an alternative geometry 
that retains the high level of blurring would be to introduce an 
outflowing vertically extended corona, where the continuum 
emission comes from large height regions and the reflection 
spectrum is produced by the emission from the base of the 
corona. However this is a complex geometry to model, mak¬ 
ing it hard to implement and beyond the scope of this work. In 
any case, the requirements for high spin and small inner disk 
radius would have to be retained in any such model, otherwise 
the broad line could not be produced. 

The weak narrow line detected in the NuSTAR and XIS1 
spectra has been_detected previously in the hard state with 
Chandra ( [Miller et al.|2002|), b ut was not detected in the soft 
state by either Chandra ( [Feng, Tennant & Z hang 2003) or 
in the combined Nu STAR and Suzaku soft state observation 
(Tomsick et al.|2014|). This difference is most likely due to the 
lower flux level in the hard state and the smaller contribution 
from relativistic reflection, which is typically much stronger 
in the soft state. The equivalent width of th e line w e fin d, 
10 d= 1 eV, is lower than that found by |Miller et al.| ( |2002| ), 
16^2 eV, but not by a great deal. This could indicate changes 
in the outer disk or stellar wind, where some of the line flux 
is thought to originate, or could be due to systematic effects 
introduced by the different levels of absorption in the observa¬ 
tions - the Chandra spectrum was taken when there was very 
little absorption by the stellar wind. 

Simple phenomenological cut-off power law models for the 
continuum spectrum of Cyg X-l fail to properly reproduce the 
sharp high energy cut-off observed. This is a known problem 
with these models (e.g. [Zdziarski et al.[2003| , as the turnover 
predicted by Comtonization can be significantly steeped than 
an exponential. We find that a simple disk black body, ther¬ 
mal Comptonization and reflection model cannot adequately 
describe the broad-band spectrum of the source. There are 
several potential explanations for this, which cannot be easily 
differentiated. The Comptonization model used may not be 
suitable for fitting the spectrum (COMPTT is only an analytical 
approximation of thermal Comptonization; |Titarchuk~||1994[ ). 
We find a significantly improved fit using the more com¬ 
plex EQPAIR model, which includes both thermal and non- 


thermal Comptonization. Alternatively, the additional curva¬ 
ture may be due to deficiencies in the reflection modelling. 
The Relxilllp model used here is calculated for a density 
of 10 15 cm -3 , which is closer to the expected density of an 
AGN accretion disk than that of an XRB. Differences in the 
density of the reflecting medium are unlikely to make a differ¬ 
ence to the parameters returned for the relativistic blurring of 
the iron line, but can cause subtle changes to the broad-band 
reflection spectrum (Fabian et al., in prep). 

We find that the same best-fit eqpair model gives a sig¬ 
nificantly worse fit to the data when we exclude non-thermal 
electrons. This, combined with the failure of the COMPTT 
model, demonstrates that thermal Comptonization alone can¬ 
not describe the continuum of Cyg X-l in the hard state. This 
leaves two possibilities - a combination of thermal and non- 
thermal Comptonizationasji^ourbest^niodehorajetbased 
model ( [Markoff, Nowak & Wilms |2005HMaitra et al.|2009[ ). 
Although we do not test jet descriptions here, we demonstrate 
that the high quality broad band observations of XRBs now 
possible with NuSTAR can distinguish between different con¬ 
tinuum models. This is a promising avenue for further re¬ 
search. If the Comptonization interpretation is correct, this 
may provide a possible explanation for the difference in the 
measured continuum spectrum and that inferred from the re¬ 
flection spectrum. The reflected emission could potentially be 
produced primarily from one population only, leading to the 
observed difference. 


5. CONCLUSIONS 

The NuSTAR and Suzaku spectra of the hard state in Cyg X- 
1 have resulted in a broad band spectrum of exceptional qual¬ 
ity, which enables weak reflection features to be measured 
more precisely than ever before. Our main conclusions are 
summarised below: 

• The broad-band spectrum of Cyg X-l shows a very 
steep cut-off above ^ 100 keV, detected with the 
GSO, corresponding to a plasma temperature of 43.4 db 
0.8 keV. This cut-off is not apparent in the NuSTAR 
spectrum alone and is worth considering in other NuS¬ 
TAR spectra of hard-state binaries, where simultaneous 
GSO data may not be available to measure the cut-off 
energy. 

• Fitting the iron line energy band alone (from 4-10 keV) 
we find that the line cannot be modelled with a combi¬ 
nation of narrow lines and ionized absorption, and that 
a relativistic line profile is needed. In addition, we find 
a weak narrow line is needed to fully describe the pro¬ 
file. 

• When we carefully model the relativistic blurring of the 
iron line, we find that it requires a high degree of smear¬ 
ing, with a high spin (a > 0.97) and small inner ra¬ 
dius (ri n < 1.7 nscoX as we ll as a low source height 
(h < 1.56 nsco)- This rules out significant truncation 
of the disk at the transition to the hard state. While 
the disk remains close to the ISCO in the bright part of 
the hard state, disk truncation at much lower Eddington 
fractions remains possible. 

• A simple comparison of the line profiles from the soft 
and hard states shows that the line has not become sig¬ 
nificantly narrower. While the high energy side of the 
line appears to drop earlier in the hard state, the red 
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wing extends to similarly low energies, and the differ¬ 
ences can be attributed to differences in the ionization 
of the disk. 

• We confirm the high level of blurring with fits to the 
full spectrum, simultaneously modelling the continuum 
and reflection spectra. After running extensive MCMC 
chains to explore the parameter space, we do not find 
any solutions with r m > 3 vq at the 3cr confidence level. 

• We find that thermal Comptonization models (COMPTT 
and EQPAIR without non-thermal electrons) cannot ade¬ 
quately describe the continuum. Instead a hybrid model 
is required for Comptonization to fit the data. 
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